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8 : Abstract 

We study gravitationally bound static and spherically symmetric configurations of k-essence 

| — > : 

fields. In particular, we investigate whether these configurations can reproduce the properties 

<N . 

of dark matter haloes. The classes of Lagrangians we consider lead to non-isotropic fluids with 

<N 

J> ■ barotropic and polytropic equations of state. The latter include microscopic realizations of the 

C-* . 

often-considered Chaplygin gases, which we find can cluster into dark matter halo-like objects 

(N 

with flat rotation curves, while exhibiting a dark energy-like negative pressure on cosmological 

in 
o 

& • initial value formulation of non-canonical scalar field theories, and we also discuss a new class of 

q de Sitter solutions with spacelike field gradients. 

in 
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scales. We complement our studies with a series of formal general results about the stability and 
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I. INTRODUCTION 



The recent remarkable progress in cataloging the demographics of our universe is leading 
us inexorably towards the conclusion that it consists of about 70% dark energy, 25% dark 
matter and 5% everything else such as radiation and regular baryonic matter The 
dominant dark energy component, which has a negative pressure, drives the acceleration 
of the universe ; while the pressureless dark matter gravitationally collapses into non- 
linear objects, forming the backbone of the large-scale structure of the luminous universe 
we see |4|. 

While the standard model of particle physics provides us a with a quite successful de- 
scription of this 5% of baryonic matter in the universe, our knowledge of both dark energy 
and dark matter is tenuous at best. One of the more promising incarnations of dark energy 
involves non-canonical scalar fields, also known as k- fields [5j]. They were first introduced 

n 

in the context of inflationary models [f|, but it was soon realized that they offer great flex- 
ibility in designing various cosmological fluids with negative pressure, generally known as 
"k-essence" (3, 0]. Examples of fluids that can be cast as k- fields include Chaplygin and 



3orn-infeld gases jsl 12|, tachyon matter jl3j |. phantom energy l^J and ghost condensates 

Conventionally, we think of dark matter and dark energy as two separate entities. How- 
ever there is a certain elegance in finding a unified model that describes both simultaneously. 
Although Chaplygin gases might embody such a unified description P, phenomenolog- 
ically successful models are hardly distinguishable from conventional ACDM, as far as the 
evolution of the background and its linear perturbations are concerned Recently, k- 

essence fluids have also found explicit application as possible ways to unify dark matter and 
dark energy in a model proposed by Scherrer |ig| . The formation of structure in Schrerrer's 
unified description has been studied by Giannakis and Hu j^J. 

Irrespective of the model under consideration, dark matter perturbations have to collapse 
to form approximately static objects which serve as gravitational potential wells where lu- 
minous matter such as stars and galaxies cluster. The standard paradigm to describe such 
large-scale structure is the halo model 21 1, which posits that cold dark matter agglomer- 
ates into roughly spherical objects called "haloes" . In this work, we extend the analysis of 



structure formation in k-essence models beyond the perturbative regime 



in partic- 
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ular we are interested in the structure of haloes made of k-essence fluids. Specifically, we 
study static and spherically symmetric solutions of the Einstein field equations sourced by a 
number of astrophysically interesting k-essence models. These include the usual barotropic 
gases considered in standard cosmology (p oc p), the polytropic gases (p oc p k , k > 0) con- 
ventionally studied in astrophysics and the generalized Chaplygin/Born-Infeld gas models 
(p oc p~ k , k > 0) proposed as forms of dark energy. The final case is especially intriguing, 
since cosmologically the fluid has negative pressure while in the static regime, as we shall see, 
it can form spherical objects with roughly flat rotation curves, implying that cosmological 
perturbations of dark energy may collapse to form dark matter haloes. 

Before we proceed, we would like to comment on some philosophical questions about our 
approach. In general, dark matter is treated as an ensemble of low momentum particles, 
as opposed to as a classical field. In other words, given a matter field Lagrangian, we 
conventionally assume that dark matter is described by an incoherent excitation of the field 
quanta, rather than by a coherent one. However, in this paper, we take the second approach: 
we directly solve the Einstein equations with k-field matter as a classical source. This is 



akin to some treatments of bosonic stars 



25] and bosonic dark matter haloes 2 



27]. 



This paper is organized as follows. In Section |Hj we state our model and derive the weak 
field equations that spherically symmetric and static configurations of a k-field have to sat- 
isfy. In Section ITTT1 we solve these equations for two particular forms of a k-field Lagrangian: 
a barotropic model (Section IlII AJ) and a polytropic/Chaplygin gas model (Section IlIIBft . 
Finally in Section IIVt we present our conclusions. In addition, we have included several 
appendices which contain various derivations of relations and asides, that would have clut- 
tered up the main text. However, we would like to highlight Appendix El where we provide 
a careful examination of the conditions one must impose upon the k-field Lagrangian for it 
to have a valid initial value formulation, and Appendix IU1 where we discuss a new class of 
de Sitter solutions with spacelike gradients (as opposed to the timelike de Sitter solutions 



of k-inflation 



)• 



II. PROPERTIES OF INHOMOGENEOUS K-ESSENCE FLUIDS 

In this section, we introduce the Lagrangian of our model. We write down the field 
equation of motion and compute its stress-energy tensor, identifying the relevant physical 
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variables along the way. For later use, we comment on the important differences between 
time-dependent and homogeneous — "cosmo logical" — field configurations and its static in- 
homogeneous counterparts. As we shall see, while the stress tensor of a cosmological field 
does describe a perfect fluid, the stress tensor of a static field does not. Following this 
discussion, we focus on the equations of motion of static and spherically symmetric field 
configurations by deriving the generalized hydrostatic equations and their weak field limit 
forms. We postpone the investigation of their solutions to the next section. 



A. The Model 

Our starting point is a non-canonical scalar field minimally coupled to gravity 

R 



S = Id xa/— g 



+ L(X) 



(1) 



16ttG 

where the Lagrangian L is some arbitrary function of the squared gradient of the scalar field 

X = ^V„¥>W- (2) 

Following 0, we call such a scalar field a "k- field". Note that we are working in the 
(+, — , — , — ) metric signature convention. In general L can be a function of both X and ip, 
though here we assume that the Lagrangian only depends on X. This purely X-dependent 
Lagrangian is the most general one compatible with a shift symmetry ip — > ip + const, and 
second order equations of motion. The shift symmetry implies that there is a conserved 
current dL/dX ■ V^y?, whose conservation equation is simply the field equation of motion 

V,. = 0. (3) 

Varying the action with respect to the metric g^ u yields the k-field energy momentum tensor 

Tfxv = ^ V ^V„^ - L{X) 9lxv . (4) 

The behavior of such a k-essence fluid depends on the scalar field gradient, which in turn 
depends on the symmetries of the spacetime the field is living in. In the often considered 
case of a k-field in a Friedman-Robertson- Walker (FRW) universe, the symmetry of the 
Robert son- Walker metric demands that the scalar field be spatially homogeneous (up to 
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Weak Energy Condition (WEC) 


p>0,p + p>0 


Null Energy Condition (NEC) 


p + p>0 


Dominant Energy Condition (DEC) 


\p\ > \p\ and WEC 


Null Dominant Energy Condition (NDEC) 


\P\ > \P\ 



TABLE I: Energy conditions for a spacelike scalar field gradient. 

small perturbations), so its gradient is timelike, X > 0. In this case, the energy momentum 
tensor has perfect fluid form, 



Tfw = {pc + P^UpUt 



Pc 9/J.U 1 



(5) 



where the energy density p c and the pressure p c are given by 

p c = L, Pc = 2X^-p c (X>0). 



(6) 



We have added a subscript c for the "cosmological" density and pressure terms to distinguish 
them from their counterparts in the inhomogeneous case below. Note that the four-velocity 
of the fluid is proportional to the scalar field gradient. In contrast, if the scalar field is 
static but spatially inhomogeneous, then its gradient is spacelike, X < 0, and the energy 
momentum is not of prefect fluid form. Instead, it can be cast as (note the last term) 



T^v = {p + p) n^n v + pg^, 



(7) 



where n u is a spacelike unit vector in the direction of the scalar gradient, n^ 1 = V^V/V - 2X. 
Hence, an observer whose four-velocity is perpendicular to n M measures an energy density p, 
a pressure in the direction perpendicular to n M equal to — p, and a pressure in the direction 
parallel to n M equal to p, where the energy density and radial pressure are given by 



p = -L, p = 2X^-p (X<0). (8) 

In other words, the pressure of the fluid is no longer isotropic [5]. Since dark matter haloes 
are roughly spherical and static, in this paper we specialize to spherically symmetric and 
static configurations of these k-fields. For later reference, we summarize the classical energy 
conditions for these field configurations in Table HI For a definition of the classical energy 



conditions, the reader is referred to \'2i 



29 
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B. Spherically Symmetric Static Solutions 



The most general static and spherically symmetric spacetime metric is 

ds 2 = e 2a ^dt 2 - e 2 ^dr 2 - r 2 (d6 2 + sin 2 ## 2 ), (9) 

where a and (3 are two r-dependent functions. Varying the action (JTJ) with respect to the 
metric leads to the Einstein equation G^ u = 8-tiGT^. For the metric Q, G^ u is diagonal, so 
it follows from the Einstein equation that T^ u also has to be diagonal. Using equation to 
compute T tr , we find that the first term of the energy momentum tensor (@J), (dL/dX)dtp d r p, 
has to vanish. Excluding the special case dL/dX = 0, we conclude that ip only depends 
either on t or r. Solutions where the scalar field only depends on time, p = p(t), have been 
considered in jjy]. In this paper, we restrict our attention to solutions where ip is only a 
function of the radius, i.e. 

p = <p(r) (10) 

Then, the non-vanishing terms of the energy momentum tensor are 

T\ = p = -L (12a) 

T\ ee p = 2^-X - L (12b) 
dX 

T% = T% = —L, (12c) 

Note the presence of a significant difference between the radial pressure p and the transverse 
pressure L. The next step in our derivation is to write down the relations of hydrostatic 
equilibrium for the system, which are the various components of the Einstein equation. 
However, for the applications we have in mind, it is not necessary to consider the full non- 
linear Einstein equations. Indeed, for most astrophysical situations the weak field limit 
conditions 

Gm 

< 1, (13a) 



r 

,2 



AnGpr 2 < 1 (13b) 
apply. Then, if we define 

2Fj ^ 2Gm n Gm , 

e = 1 — , or p w — — , (14) 
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we find that the Einstein equations yield the following set of structure equations 



dm 
dr 
da 

dr 
dp 
dr 



47T pr 2 , 

Gm + 4irGpr 3 



r(r — 2Gm) 



[2 da\ . , 2. . 
{- + Tr ) ( P + p)*- r (p + p), 



Gm + AnGp r 3 



(15b) 



(15a) 



(15c) 



where the equal sign denotes the exact relations, while the approximation denotes their 
weak field regime cousins. The first two equations are the tt and rr components of the 
Einstein equation respectively, while the 3rd equation is a combination of the rr and the 99 
components. These "master" equations are very similar to the Tolman-Oppenheimer-Volkov 
equations of stellar structure |32|, where matter is assumed to be a perfect fluid. Since the 
k-essence fluid is not perfect, there is an additional 2/r term in equation ()15c|) arising from 
the tangential pressure terms. Note that (|15c|) is equivalent to the equation of motion (J3J). 
Thus, this set of equations suffices to determine all the unknowns in the problem, which are 
a, (3 and p = — L (p is related to L in terms of equation (JSJ)). 

III. K-ESSENCE FLUIDS AS DARK MATTER HALOES 

Since the Lagrangian L(X) is an arbitrary function, the space of possible field configura- 
tions which satisfy the Einstein equations is infinite. The only statement we can make about 
the fluid at this moment is that it is "isentropic" , by which we mean that the pressure is 
simply a function of its density p = p(p). We are, however, interested in functions that can 
lead to dark matter halo-like configurations. To this end, we focus on two particular forms 
of the Lagrangian which are both simple enough to solve and have such properties. As we 
discuss next, our choices of Lagrangians can be also motivated as follows. 

The macroscopic behavior of a fluid is described by a set of state parameters. This 
set of parameters determines its stress energy tensor, which is then coupled to gravity via 
Einstein's equations. These parameters are in general not independent, but constrained 
instead by an equation of state, which depends on the microscopic properties of the fluid. 
In many cases, the equation of state can be written in the following form 



p = Kp 1 



(16) 
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where K and 7 are constants. If 7 = 1, we call such a fluid barotropic, and if 7 7^ 1 we call it 
polytropic. An example of an isotropic barotopic fluid is a bath of photons where K = 1/3, 
while an example of a polytropic fluid is a bath of non- relativist ic degenerate electrons of a 
white dwarf star, where 7 = 5/3. Both the above examples, and other isotropic fluids, have 
been extensively studied in the field of stellar structure and evolution (3^]. In this paper, we 
investigate the structure of non-isotropic k-essence fluid haloes. In analogy with equation 
(fTfiJl we postulate the equation of state for the k-essence fluid to be of the form (fTBj) where 
p and p are the pressure and energy density defined in equation (jHJ). Then, it follows that 
for homogeneous configurations of the scalar field, the equation of state in a FRW universe 
(if defined) has the inverse form, 

[P « Plstatic ~ [Pc OC P\ h ] homogcneous • (17) 

One can see this simply by comparing the relations between p (p c ) and p (p c ) in the two cases, 
namely equations (jOJ and (jHJ), and noting that the roles of p (p c ) and p (p c ) are interchanged 
between the two. Furthermore, we show in Appendix El that this inversion also applies to 
the squared sound speed of its perturbations c 2 s = Sp/5p, as one would naively expect. 

Given the form of the equation of state (|16p. we can work backwards and derive the 
k-field Lagrangian we need to describe such a fluid. It turns out that the corresponding 
Lagrangians span simple polynomials or non-analytic functions already proposed in quite 
different contexts. We show in Appendix that, regardless of the Lagrangian, the only 
solution to the equations of motion regular at the origin is de Sitter space. Therefore, the 
solutions we discuss next are not well-defined for all values of r, and they generically become 
singular as r — > 0. Since we expect new physics as this limit is approached, we shall not be 
concerned about this fact. In any case, because we only study the weak field limit of the 
gravitational equations, we do not encounter this singular behavior in the regime of validity 
of our approximations. 

The solutions we find this way are characterized by two boundary conditions, the energy 
density po an d the mass m at a given, arbitrary value of the radius r . Presumably, their 
value is determined by the environment that lead to the formation of the halo, a process 
we cannot address with our static solutions. Since we cannot determine the value of these 
parameters, and because they determine some of the halo observables, like their mass, we 
mainly focus on the scaling of our solutions and on their main properties. 



S 



A. Barotropic Lagrangians 



We first consider a barotropic k-essence fluid p oc p. Suppose that the Lagrangian L is 
an analytic function of X. Then, L can be expanded in a power series, 

L = Ml T,^(wT' (18) 

where the c n are dimensionless coefficients and M is an arbitrary mass scale. For a given 
value of X, we expect one particular term of the series to dominate, depending on the values 
of c n . Hence, we can assume that (at least for a certain range of values of X) the Lagrangian 
is given by the monomial 

l = m'(wJ- < 19 > 

If the scalar field is spatially homogeneous, as it must be if it shares the symmetries of a 
FRW spacetime, such a scalar field behaves as a perfect fluid with a barotropic equation of 
state 

p c = wp c , (20) 
where for convenience we have defined the equation of state parameter 

In contrast, as we have described in Section |Hj for static and spherically symmetric config- 
urations the inverse relation holds 

p = £ (22) 
w 

Similarly, inserting the Lagrangian (fT9~j) into the sound speed equation (jB8|) . we find that in 
the static case the sound speed is given by 

c 2 s = - (23) 
w 

as expected from the equation of state (|22[). The speed of sound is superluminal if w < 1 
(n > 1.) As detailed in Appendix El this is not a violation of causality, since the k-field 
propagates on the "light-cone" of a different metric. In fact, unless w — 1, the speed is 
superluminal for either homogeneous or static field configurations. Equation (}2*3*j) also shows 
that if w < (n < 1/2), the k-field is unstable to perturbations. Notice that in order for 
the energy density to be positive in both cosmological and static cases n has to be odd. 



In order to compute the energy density profile of a halo, we insert equation (}2*2*j) into the 
weak field structure equations (JT3j), and find that the energy density is a simple power law 



where po and mo are the energy density and mass at tq. The energy density p is posi- 
tive if po is positive. Now we divide the solutions into three different classes: n < 1/2, 
l/2<n<3/2 and n > 3/2. We will discuss each in that order, with the last case being the 
most astrophysically interesting. 

In the first case, n < 1/2, the equation of state parameter w is negative, leading to a 
negative squared sound speed cj. < via equation (J2l!j) . Thus, configurations with w < are 
unstable, and we shall not consider them here except to mention that they do not violate 
the weak field conditions as r — > 0. 

Before we proceed to the next range, we note as an aside that w = oo (n = 1/2) is a 
special case; using equation we see that the radial pressure p identically vanishes for 
this n. Although this is not compatible with the weak field equation (|15cjl . it is a valid 
solution of the full equations of motion if da/dr = —2/r. In that case the equations of 
motion are exactly solvable, and we find that r is a timelike variable. Hence, this solution 
does not describe a static configuration. 

The second case, 1/2 < n < 3/2, implies that 1/2 < w < oo, which is the only range 
of w where the mass of the halo remains finite as r —>■ oo. Note that in order for the weak 
field approximation to remain valid as r — > oo, w has to be positive (n > 1/2), which means 
that all finite mass solutions can be made consistent with the weak field approximation by 
an appropriate choice of initial conditions. On the other hand, if we try to extrapolate our 
solutions to small radii, we find that the week field limit is violated as r — ► 0. 

Finally, we discuss the case where n > 3/2, implying that w < 1/2. This class of solutions, 
similar to the 1/2 < n < 3/2 case above, violates the weak field limit as r — > 0. The mass 
of their haloes also diverges 1 as r — > oo. As described above, the equation of state of our 

1 This does not mean that they are unphysical however. As we mentioned above, if the Lagrangian is a 



m r, 




(24) 



(25) 
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barotropic fluid in a cosmological space time is p c = wp c . Thus if n is very large, w is close 
to zero, which is simply the equation of state for cosmological dark matter. Furthermore, 
the cosmological sound speed for this matter is c 2 s = w which also approaches as n tends 
to large values. Hence, in the linear regime this Lagrangian exhibits all the properties of 
plain vanilla dark matter. In a dark-matter dominated universe, linear perturbations grow 
with the scale factor. Once they reach the non-linear regime they collapse into dark matter 
haloes, which we identify here with the static spherically symmetric matter configurations 
we have found. As we have shown in equation ()24|) . for large n their energy density profile 
at large r becomes 

P ~ Po (yj (26) 
which is just that for a singular isothermal sphere. The singular isothermal sphere is a good 
approximate model for a galactic halo, because its rotation curves 2 flatten out as the radius 
tends to infinity, as consistently observed in galaxy rotation curves js^J. Since the observed 
luminous mass is insufficient to explain such a velocity profile, this fact is usually attributed 
to the existence of non-luminous (cold) dark matter. 

Rotation curve observations are usually measurements of the Doppler shift of the light 
emitted by HII regions or, for nearby galaxies, stars. As shown in Appendix the apparent 
rotation velocity of an object at radius r inside the k-essence halo, v (r), is given by equation 

dun, 

v 2 = Gm(r) +47rG7>2 
r 

For the barotropic model, the right hand side of equation (j2*7|) is analytically solvable. By 
plugging in the density profile equation into the equation above, we find 

v\r)*4irGr 2 p Q ^-^-(^) , (28) 
w — 2w z \r J 

which as w — > (n — > oo) becomes a constant, i.e. the rotation curve becomes flat at large 
r. While this is an interesting model of dark matter, as it exhibits the dark matter-like 



polynomial in X, rather than a monomial, different monomials dominate at different values of X. At 
large values of X, high powers are relevant, and for low values of X, low powers are dominant. Now, 
because the energy density monotonically decreases as the radius r increases, far from the halo center 
we expect the Lagrangian to be dominated by a low-power monomial (typically the linear term, n = 1.) 
These low-power terms then render the mass of the halo finite. In other words, as r — > oo, one can expect 
the solution to return to that of the second case. 
2 i.e. the circular velocity of objects, say stars, orbiting around the galactic center. 
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n < 1/2 


1/2 < n < 3/2 


n > 3/2 




w < 


1/2 < w < oo 


< io < 1/2 


mass 




finite 


diverges 


valid weak field regime 


r -> 


r — » oo 


r — > oo 


c 2 


< 


> 


> 



TABLE II: Properties of barotropic solutions. 



properties cosmologically and in the non-linear regime, it is nevertheless not a very natural 
model, since n has to be extremely large. Indeed, if the sound of speed dark matter is non- 
vanishing, density perturbations are prevented from growing by the non- vanishing pressure, 
which causes density perturbations to oscillate instead. The absence of these oscillations in 
the (galaxy) power spectrum down to scales of around 1 Mpc j^] hence restricts the speed 
of sound of dark matter to be less than about 10~ 3 , which translates into n > 10 6 . 

The properties of barotropic solutions are summarized in Table ITT1 Note that the solutions 
go into the strong field regime either as r — > or r — ► oo. 



B. Polytropic Lagrangians 

In this subsection we consider a class of non-analytic Lagrangians that lead to a polytropic 
equation of state 

7 

(29) 

where 7, and p* > are constants. By inserting equation (ffilj) ii 
finds that such k-essence fluids are given by a Lagrangian of the form 

1 

1-7 

(30) 

The first minus sign (in front of the bracket) guarantees that the halo energy density is 
positive, while the second one (inside the parenthesis) guarantees that the Lagrangian is 
well-defined for spacelike X and arbitrary 7. From this choice of signs it also follows that 
L x is positive whenever the expression inside the brackets is positive. As shown in Appendix 
iBl this is one of the conditions for the quantum stability of the configuration. If in addition, 

12 




;o equation 




one 



-M 4 X)^ - % 
pi 



Cg > 0, the configuration is both classically and quantum- mechanically stable. In order to 
compute the speed of sound, we insert the Lagrangian (j3T?j) into the sound speed equation 
(|B8|) and find that the sound speed in the static case is 

cl = l V - (X<0). (31) 

Since by assumption and construction p > 0, this means that the sign of jp must be positive 
in order for the perturbations to be stable. This result allows us to rule out many of the 
possible cases in our following discussion of polytropic solutions. 

The class of Lagrangians (j3*Uj) is interesting for several different reasons. For instance, 

3811 . The Born-Infeld 



if 7 = — 1 and < 0, it becomes the scalar Born-Infeld Lagrangian 
Lagrangian has found many applications in high ener gy physics and cosmology |3jJ]. It 
naturally appears in the description of moving branes (ill ] , a nd in particular, it has been 
argued to describe the dynamics of open string tachyons 
suggested to behave as a pressureless fluid, i.e. as dark matter 



|| , a i 

n 
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he latter have been also 
If the scalar field is 



homogeneous, the Lagrangian (j3~Uj) leads to the polytropic equation of state 

Pc = (-/>•) (—\ H < 0. (32) 



-P*. 

Insisting that the energy density be positive leads to a negative cosmological pressure, a key 
requirement for an accelerating universe. Note that the speed of sound in this case is 



c\ = 1^ (X > 0), (33) 

which is dark matter- like (small and positive) if 7 <C — 1, even when p c « — p a . 

For negative 7 (which includes the Born-Infeld case), the equation of state of a homo- 
geneous k-field (|32|) agrees with the one of a generalized Chaplygin gas. Chaplygin gases 
have been widely considered as a unified description of dark matter and dark energy |9j. An 
advantage of our approach is that the Lagrangian (J3U)) provides a microscopic description 
of a Chaplygin gas. In particular, it shows that in some cases is not sufficient to specify 
the equation of state of a homogeneous component; one must specify the full Lagrangian 
to obtain a description of all its possible configurations. In our realization, the Chaplygin 
"gas" does not behave as a perfect fluid when the k-field is static. 

Inserting the equation of state (J2H|) into the equation of motion ()15c|) we find 



P = Po 



2-27 

Po \ I r \ ~< po 



1 + 

Po) \r J Po 



7-1 

(34) 
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where po and po are the pressure and density respectively at r . Note that the signs of po 
and po must follow that of p* and p* respectively The expression inside the parenthesis of 
equation (J53|) vanishes at the radius 



The behavior of the density profile depends on the the position of this point, which depends 
on the values of po and 7; say, if p(rg) vanishes, then r$ is a "boundary" of the halo. Since 
Po fixes the overall sign of the pressure, equation (JHTJ, stability requires that 7P0 > 0. On 
the other hand, the sign of p is determined by p*, so that the Lagrangian itself (as opposed 
to the initial conditions) is what fixes the properties of the solutions. In the following we 
discuss the different relevant cases, {p* > 0,7 < 0} and {p* < 0,7 > 0}, separately. The 
results are summarized in Table IIII1 

a. p* > and < 7 < 1 The energy density has a pole at < r , and as r — ► 00 
it decays as p oc r~ 2 / 7 . Hence, the total mass of the halo is finite for 7 < 2/3. Note that 
the ratio p/p — > 00 as r — > 0, violating the energy conditions listed in Table HI However, 
these solutions are both classically and quantum-mechanically stable, so they are perfectly 
sensible scalar field configurations. The Lagrangian for this case is plotted in Fig. ^ while 
its density profile is plotted in Fig. |21 



FIG. 1: A plot of the Lagrangian for > and < 7 < 1 (case a). The zero occurs at 




(35) 



L[X] 




M 4 X 



X = —M 4 {p * I ' pl) 2 ^ 1 J \ and the pole of p{rs) corresponds to the X 



00 limit. Note that 



the Lagrangian is not defined for positive values of X. 
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FIG. 2: The density profile p for a model with > and < 7 < 1 (case a). The parameters are 
p*/p* = 0.9 and 7 = 1/2. Since 7 < 2/3, this halo has a finite mass. 

Again, we test such models against observations of the motion of objects around a galactic 
halo by computing the rotation curves of test particles. At large r the density of such a 
model becomes 

/ n \ V(7-l) / r \ -2/7 

© • (36) 

Thus, in order to calculate the behaviour of rotation curves at large r, we simply insert 
equation (j3T?j) into equation (jA15|) . 



(37) 



Focusing on the first term, we see that the mass of the halo grows at most linearly in r, so 
this term is constant or decays with r. Meanwhile the second term of equation (|37jl is 

/ „ \ 7/(7-1) 

AnGr 2 p = AnGrlpo ( 1 + f (38) 
which is a constant. Hence, the rotation curve v(r) approaches 



r V Po/ 



7/(7-1) 



(39) 



which is constant for sufficiently big r. As mentioned above, the flatness of stellar rotation 
curves is conventionally assumed to be due to the presence of dark matter which, by def- 
inition, is pressureless. In our halo, the contribution of the energy density to the rotation 
curve vanishes or becomes a constant at large r, and it is the non- vanishing radial pressure 
what keeps the rotation curve flat. 
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b. p* > and 7 > 1. The Lagrangian and density profile corresponding to this case are 

shown in Figs 01 and 0] respectively. For these values of the parameters, the energy density 

vanishes at r$- At this point of vanishing energy density, the solution can be joined to the 

vacuum, ps = Ps = 0. Therefore, this halo has a surface, and it seems more appropriate to 

call it a "k-star" . Spherically symmetric configurations of (canonical) classical scalar fields 

have been widely considered in the literature under the name of "boson stars" (see e.g. 2sj|-) 

The difference between these solutions and ours are that in the former the energy density 

smoothly approaches zero at infinity, while in the latter the energy density vanishes at a 

finite radius. We shall not analyze the properties of these stars here, since the only relevant 

observable, their mass, depends on the initial conditions, especially the core mass. However, 

these configurations could still pla y a secondary role in the dark matter puzzle as massive 

PI 

compact halo objects (MACHOS) |4l|. 



l[x; 



M 4 X 



FIG. 3: Figure showing the Lagrangian (j30j) for > and 7 > 1 (case b). The Lagrangian has a 
pole at X = —M~ i (p*/p2) 2 /( 1 ~~ / \ and a zero at X = 0. Note that the Lagrangian is defined for 
positive values of X only if 7 is an odd integer, with the resulting plot being a mirror image of 
X < 0. 



c. < and 7 < 0. Finally, we consider the case of a generalized Chaplygin gas, 
which has been suggested as a unified description of dark matter and dark energy. Since the 
Lagrangian (PU|) implies that po/po > — 1, the energy density has a pole at rs < r , and as 
r — > 00, it approaches the constant (see Fig. 

p2 



P* 



(40) 
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FIG. 4: Figure showing the density profile for a model with > and 7 > 1 (case b). The sample 
parameters are p*/p* = 3/4 and 7 = 3/2. There is a point of zero density at finite r, and the 
solution can be joined to the vacuum at that point. 

Therefore, in this form solutions asymptotically approach de Sitter spacetime (one can also 
check that the pressure approaches —A), which is what we expect, since Chaplygin gases 
lead to cosmic acceleration. From the point of view of the Lagrangian, which we plot in 
Fig. this is easy to understand. As shown in the Figure, at X = the Lagrangian is 
negative and linear in X, and hence, it is not surprising to find solutions that asymptotically 
approach de Sitter. In order to disentangle the vacuum energy from the halo density, we 
substract this cosmological constant from the Lagrangian, L —>■ L + A , where L is our 
original polytropic Lagrangian. Then, the total density and pressure become p = p — A and 
p = p + A respectively, where tilded quantities are associated with the original Lagrangian L. 
Because this shift leaves the weak field equation of motion (j!5c|) invariant, the new density 
profile p is simply the original profile p shifted by the same constant, 



Thus one can see that by adding the constant A to the Lagrangian, we have renormalized 
the asymptotic value of the energy density at infinity to zero. 

In order to compute rotation curves for this halo, we expand the root around its asymp- 



P 



Pa 




(41) 
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M 4 X 



FIG. 5: Figure showing the Lagrangian ()3U|) for p* < and 7 < (case c). We plot both the 
original Lagrangian (full line) and the renormalized Lagrangian (dashed line). The asymptotic 
value of p{r) corresponds to the energy density at X = 0, which is a non-zero constant (full line). 
Note that the Lagrangian is defined for positive values of X only if 7 is a negative odd number. If 
(1 — 7)/2 is even, the plot of the Lagrangian for X > is a mirror image of the one for X < 0. 
On the other hand, if (7 — l)/2 is odd, the Lagrangian is differentiable at X = and vanishes at 
a finite positive value of X. 
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FIG. 6: Figure showing the density profile for a model with < and 7 < (case c). The 
sample parameters are p*//?* = —0.4 and 7 = —1/2. The density profile asymptopes at large r 
to a constant (full line), which can be subtracted off by adding a cosmological constant to the 
Lagrangian (dashed line), as we do in the text. Note that the r/ro axis begins at unit value. 
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totic value, obtaining 



P = Po 



1 + ^ 



Po 



i 

7-1 



1 - 7 V Po J \ -Po 
Similarly, the pressure term is approximately given by 

Po 



P = Po 



JL_ 1 + Po 

1 - 7 V Po 



-Po 



i 

7-1 



2-27 

1 



> 0. 



(42) 



2-27 



> 0. 



(43) 



Note that the total pressure is positive, though p is negative. Recall now that the rotation 
curve is given by equation (jA15|) 



= Gm ( r ) + 47r £r 2 p(r). 



(44) 



The pressure scales like r 2 p ~ r 2//7 , and the mass is m(r) ~ m + mir 1+2//7 , where uiq and 
mi are integration constants. Therefore, the rotation curve has the form 

Gm 



v (r) 



+ iVr 2/7 , 



(45) 



where iV is a constant too. Hence, if 7 > —2, the rotation curve goes as 1/r, and if 7 < —2 
the rotation curve decays as r 2 ^ at large r. 

It follows from our considerations that polytropes with 7 — 1 are particularly inter- 
esting. On one hand, their rotation curves are flat, so these models correctly reproduce 
the structure of dark matter haloes. On the other hand, their cosmological sound speed is 
also dark matter-like, and it turns out that unified models based on these Chaplygin gases 
successfully reproduce the linear structures we observe in the universe 





p* > 


p* < 




< 7 < 1 


7 > 1 


7 < 


rs 


< r 


> TO 


> r 


P(rs) 


pole 


zero 


pole 


r — > 00 




p -» f- 2 /T 


constant 



TABLE III: Polytropic solutions. 
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IV. CONCLUSIONS 



We have investigated static and spherically symmetric solutions to the Einstein field equa- 
tions coupled to k-essence. Our results shed light on the formation of non-linear structures 
in models where matter can be microscopically described by a non-canonical scalar. In con- 
trast to homogeneous k-essence fluids, the energy momentum for these configurations does 
not have perfect fluid form. Rather, it describes an anisotropic fluid with different pressures 
in the radial and tangential directions. Therefore, the propagation speed of perturbations 
around these solutions differs from their cosmo logical counterpart. We have computed the 
sound speed of the former and formulated the necessary conditions for these solutions to be 
both classically and quantum-mechanically stable. 

By making simple assumptions about the form of the k-essence equation of state, we have 
found solutions to several models of k-fields which exhibit dark matter halo-like properties, 
namely, flat rotation curves. The corresponding Lagrangians are simple (as in the barotropic 
case), and physically motivated (as in the Chaplygin gas/polytropic case). In both cases, 
rotation curves are flat when the cosmological speed of sound of k-essence tends to zero. The 
latter include proposed unifications of dark matter and dark energy based on Born-Infeld 
and generalized Chaplygin gases. Among these, the ones that yield flat rotation curves are 
precisely the ones consistent with the way linear structures forms in our universe 

Our solutions are valid only in the limit where the gravitational field of the halo is weak; 
as one approaches the center of the halo, the weak field limit is generically violated as the 
density diverges. In other words, haloes made of this form of k-essence fluid also contain 
a central "cusp", whose existence in real haloes is still under debate (see e.g. 42j.) In 
any case, this is a generic problem in common cold dark matter models, and is conceivably 
resolved only when we add baryonic matter to our model. 

Dark matter is a specter in search of an identity. Although simple and appealing, the 
assumption that dark matter is made out of self-interacting particles requires involved nu- 
merical simulations to extract predictions about non-linear structure formation. In this 
paper we have explored an alternative — whether dark matter could be described by a non- 
canonical scalar field. As we have shown, this explanation has its own advantages. It makes 
simple, direct and straightforward predictions that can be immediately compared with ob- 
servations. 
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APPENDIX A: ROTATION CURVES 

Rotation curves of dark matter haloes are determined by measuring the Doppler shift of 
light emitted by the orbiting objects, such as stars or gaseous ("HII") regions. In order to 
compare our haloes with these observations, we study how the light emitted by these objects 
at radius r is red or blue-shifted. In this appendix, our goal is to derive the rotation curves 
of massive particles in circular orbits around k-essence haloes (see also 

Consider the motion of massive test particle, say a star, in a such a halo. Its trajectory 
is then described by a curve x M = (t, r, 6, <p) parameterized by some affine parameter; here 
we use its proper time r. Its four velocity is then simply = dx^/dr. Due to spherical 
symmetry, we can assume without loss of generality that the star's ecliptic is located in the 
6 = 7r/2 plane. Since the star is a massive particle, its norm is = 1, which becomes 
the constraint equation 

e 2a^2_ e 2^2_ r y = lj (A1) 

where a dot denotes a derivative with respect to proper time r. Since the metric does not 
explicitly depend on 6, the star's angular momentum J is conserved, 

J = r 2 0. (A2) 

Similarly, the metric does not explicitly depend on t, and there is a conserved energy E, 

E = e 2a i. (A3) 

Substituting equations (jA3J) and (jA2j) into equation (jAlj) . we find a first integral for the 
motion of the star, 

If 2 + V(r) = 0, (A4) 

where its effective potential is 

V(r) = (l + " \E 2 e~^\ (A5) 
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Note that the potential explicitly depends on the energy. Stationary orbits at radius r e exist 
if V and dV/ dr vanish at that radius. The former condition yields 

1 + 4 = £V 2a(re) , (A6) 

whereas the latter gives us 

J§. U + _ * + E > ( d " + f] e -*M = o. (A7) 
dr \ r e J r e V " r dr J 

Substituting equation ()A6|) into equation ()A7|) and using the 2nd TOV equation (JT3J), we 
get the following equation 

which directly relates the angular momentum J e to the density profile of the halo. 

With the angular momentum in hand, we turn to the task of finding the red or blue- 
shifting of light emitted by an object a such an orbit. Consider a photon with momentum 
p^. emitted by an orbiting star at radius r e , with its three-momentum assumed to be parallel 
to the star's three-velocity, 

/ e ,0,0,e a(re) ^ 
r e 

Since the star can move towards us or away from us, its four-velocity is given by 



rt=[£,0,0,e'* r ' ) f). (A9) 



i//< 



^e" 2a(re) ,0,0,±^ . (A10) 



Meanwhile, the frequency of the emitted light v e is given by v e = p^u^, and thus, by using 
equations (|A9)1 and ()A10|) we find the zero component of the photon momentum is 

pi = u e .(E T ■ (All) 



Let us now turn to the computation of the frequency measured by a stationary observer far 
away from the halo, at r^. Since she is stationary, her four- velocity is v M = (e~ a (' r °°\ 0, 0, 0), 
so this observer measures the photon frequency to be = p^v^ = e^^'p^. To determine 
p^,, we use the conservation of p t = e 2a p*, which yields the red/blue-shifted frequency the 
observer would measure 

Voo = Vee «<r.)- a <roo) (^E e -^re) zp £}j . ( A12 ) 
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This frequency change includes that of the photon climbing out the gravitational potential 
(i.e. its gravitational redshift); to extract out the change in frequency due to the motion 
of the star, we simply compute the difference in the observed frequency as the star moves 
towards and away from us. Using equation (jA6j) we find that to lowest order in small 
parameters this is 

_ oo_ _ _^ Al3 

v e 2z/ e r e 

which an observer interprets as being due to the peculiar velocity v of the star, 

v = — = -. (A14) 

Using equation (jA8|) To first order in J/r and dropping the subscript e, the rotation curve 
(jA14|) becomes 

Hence, once we have a solution for the energy density p, and a relation between p and the 
pressure p, we can simply insert them into the right hand side of the above equation to 
obtain the rotation curve. 

APPENDIX B: STABILITY AND INITIAL VALUE FORMULATION 

In the previous sections we have studied several different static configurations of spher- 
ically symmetric non-canonical scalar fields. Presumably, these are the end-result of the 
gravitational collapse of k-field matter. If so, they should be stable, since otherwise they 
would not be the end stage of a dynamical process. The stability of scalar field configu- 
rations depends on the microscopic properties of the scalar, i.e. on its Lagrangian L. In 
the context of cosmological solutions (where the scalar field gradient is timelike), classical 



stability requires that the speed of sound of the k- field be positive 



44), 



c - = £ = iL/IxT^k^ > (A ' >0) ' < B1) 

Note that this is a necessary (but not sufficient) condition, since in principle we can still 
suffer from instabilities in the long wavelength regime 3 . 



We thank Raul Abramo for bringing this point up. 
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On the other hand, at the quantum level, stability requires the Hamiltonian of field 
perturbations to be positive definite which translates into 

§>° and § + 2 - Y ^>° <*>°>- < B2 > 

If the last conditions are not satisfied, quantum fluctuations around the background cos- 
mological solution have negative energy, and the vacuum can decay into positive energy 
gravitons and negative energy scalar field perturbations 29]. Note that quantum stability 
implies classical stability. 

In the first subsection of this appendix, we analyze both classical and quantum stability 
of static, spherically symmetric configurations of a k-field. Again, we find that quantum 
stability implies classical stability. Moreover, it turns out that the stability conditions we 
find are intimately related to the existence of a well-posed initial value formulation, which 
we study in the second subsection of the appendix. As opposed to studies of stability, our 
approach here is quite general; we do not make any assumption about the nature of the 
scalar field or the spacetime metric, and we do not rely on perturbation theory. 



1. Stability 

Consider a static spherically solution X (r) of the equation of motion (jSJ). We would like 
to find out whether this solution is quantum-mechanically stable, i.e. whether it can decay 
through quantum processes. Because any arbitrary spacetime is locally flat, we can assume 
that locally the metric is Minkowski. Moreover, because we are studying static solutions the 
gradient field is spacelike, locally we can also assume that the field gradient points in a given 
spatial direction, say, the z direction, ip = \J —2X§z. Let us perturb this local expression, 

ip = ^-2X z + Sip, (B3) 

and study the resulting Lagrangian. Using the definition (J2J) we find that the perturbed 
gradient is 

X = X - ^-2X d z 5<p + -d^tpd^ip, (B4) 

where we are using the Minkowski metric. Hence, expanding the Lagrangian L(X) to 
quadratic order in 5<p we find 

5L = 1 [L x (X )S^ t - L x (X )(5<p% + <Vf,) - (2X L XX (X ) + L >X {X ))5 V 2 Z ] , (B5) 
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where a subscript , i denotes a partial derivative with respect to the coordinate variable 
i. The Hamiltonian of the scalar perturbations is positive definite only if the signs of the 
quadratic derivative terms in the Lagrangian are positive for time derivatives and negative 
for spatial derivatives. Thus, the condition for quantum mechanical stability is the same as 
in a cosmological background, equations (|B2|) . 

Although equation (jB5|) already suffices to determine the propagation speed of classical 
perturbations, it is more informative to consider an approach that maintains the symmetries 
of the background. Hence, we expand a scalar field perturbation 5ip in spherical harmonics 

6<p = y %26<p lm (t,r)Y lm {6,<f>), (B6) 

l,m 

and instead of expanding the Lagrangian to quadratic order, we linearize the equation of 
motion. Because of linearity, it is sufficient to consider one spherical mode at a time. 
Moreover, by rotational symmetry, we can always assume m = 0. Since matter is coupled to 
geometry through Einsteins equations, linear perturbations in the scalar field should source 
metric perturbations. However, we do not expect the speed at which small perturbations 
propagate to depend on the strength of such coupling, Newton's constant G. So we simplify 
our problem by considering the limit G — > 0. 4 In this limit matter perturbations do not 
source any metric perturbations, and one can simply study the propagation of the scalar 
field in a given, fixed, gravitational background. In this way, substituting ip = v?o( r ) + ^V 9 
into the equation of motion (jSJ) and keeping terms linear in Sip we arrive at 

1 d 2 5i Pl 2XL, XX + L x (PSw £ [r 2 e a ^(2XL xx + L x )\ d&pi 1(1 + 1) 
e 2a dt 2 L x e 2 Hr 2 eP[r*e<*-PL x ] eHr r 2 Vl 

(B7) 

Note that a static observer at radius r measures proper time dr = e a dt and proper dis- 
tance ds = e^dr. By inspection one can see then that the speed of propagation of radial 
perturbations is 

c 2 s = 2XL ' X T X + L ' X (X < 0). (B8) 
L )X 

This is the inverse of squared the speed of sound in a cosmological setting (where X is 
timelike), equation (jBl|) . Indeed, the definition of the squared speed of sound c 2 = dp/ dp 

4 Strictly speaking only limits of dimensionless quantities are physically meaningful. By G — > we mean 
the limit where every dimensionless quantity proportional to G goes to zero. 
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agrees with equation (|B8|) if the field is spacelike and p and p are given by equations (JHJ), and 
it is the inverse of equation (jB8|) if the field is timelike and p c and p c are given by equation 
(JBJ). The linearized equation of motion also shows that for non-radial perurbations the speed 
of sound is c 2 = 1, which, up to a sign, also agrees with the expectation c 2 = dp/ dp. Note 
that the sound speed is given by equation (jB8|) even if the Lagrangian explicitly depends on 
ip, i.e. L = L(X, if). 



2. Initial Value Formulation 



In this part of the Appendix we determine for what functions L the scalar field equation 
possesses a well-posed initial value formulation. After all, we are mainly dealing with solu- 
tions of a classical field theory, and requiring the equations to have a well-posed initial value 
formulation amounts to self-consistency of our approach. As mentioned above, we will not 
make any assumption about the nature of the scalar field or the metric, and we will not rely 
on perturbation theory. 

We begin by rewriting the k-field equation of motion Q as 

h^V^Vvip = 0, (B9) 

where the "induced" metric is 



dT d 2 T 

,r = Ix' r + I>o v ^ (B10) 



so apply if the 
) that if the 



For simplicity we assume that L only depends on X, though our results al 
Lagrangian depends on (p itself. It can be proven (see Theorem 10.1.3. of 
metric h^ v is Lorentzian, the field equation admits a well-posed initial value-formulation (at 
least locally), so all we have to show is that h^ v has a negative determinant. 

Let h, g and n denote the matrices with elements W, g^ u and V'VV'V respectively. 
We then have 

, fdL d 2 L \ fdL\* , / d 2 L/dX 2 , \ 

det h = det ^-g + — nj = {- j • det g • det ^1 + -^fT^ , (Bll) 

where g _1 denotes the inverse of g, and we have assumed dL/dX ^ 0. Using the identity 
det e A = e trA we find 



d 2 L/dX 2 , , 
detll + ^^g-nHexp 
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d 2 L/dX 2 -, 
1l,ll < l + ^X g n 



(B12) 



and expanding the logarithm in a Taylor series around the identity we obtain 

, / d 2 L/dX 2 ,\ ^ (-l) fc+1 fd 2 L/dX 2 \ k , i , k , ( v d 2 LldX 

trln ( 1 + ^W g n J=? * Ul75tJ tr(g n) ^ ln ( 1 + 2X ^ 

(B13) 

where we have used tr (g _1 n) fc = (2X) fc . Finally, substituting equation (jB13|) back into 
equation (|B12j) we arrive at our desired result 

, , fdL\ 4 , dL/dX + 2Xd 2 L/dX 2 

det h = (- j ■ det g • ^ ■ (B14) 

By assumption, the metric of spacetime is Lorentzian, so that det g < 0. Therefore, h is 
Lorentzian if and only if 

A ' = 2Xd^L/dX^ + dL/dX > °- (B15) 

So far, we have found that the metric h^ v is Lorentzian if condition (|B15J) is satisfied. 
Because the assumptions of the theorem mentioned above apply, the k-field equation of 
motion is guaranteed to have a well-posed initial value formulation, at least locally. An 
ingredient of such a well-posed formulation is the "causal" propagation of the scalar field. 
Note however that the metric that enters equation (jB9|) is not the space-time metric , 
but the derived metric h^ u . Therefore, causality in the propagation of the k-field means 
that sign* propagate inside the hght-cone of the metric h, and not „ Q. This observation 
sheds some light on a property of non-canonical fields that was noted early on |4£|]. If the 
Lagrangian L is an arbitrary function of X it is possible to obtain superluminal sound- 
speeds, c 2 > 1. More generally, k- fields are natural implementations of bimetric theories, 
and they might offer a consistent theoretical framework to study and implement some of the 
ideas about causality as suggested in 4f|. 

To further illustrate the causal structure tied to the propagation of a k-field, let us consider 
the inverse of the matrix h, whose components we denote by h^, 

r 1 f d 2 L/dX 2 „ „ \ . . 

^ = dLjdX ^ ~ dL/dX + 2Xd 2 L/dX 2V ^ V ^ ) ■ (B16) 

Note that the inverse of is not given by h^ v = g^ p g va h pa . The metric determines 
the causal structure of the solutions to equation (|B9|) . and the theorem mentioned above 
guarantees that no signal can travel outside the light cone defined by h pv . The question 
we are interested in is the following. Suppose that a vector field v M is null with respect to 
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h (h^uV^v^—O), is timelike, spacelike or also null with respect to gl If v M is null with 
respect to h^, it follows from equations (jB15|) and (|B16|) that 



g^v v 



IX 



(B17) 



Hence, there are two cases, as in Fig. 

• Timelike gradient (X > 0) 

In this case the vector v M is timelike for A 2 < 1 and spacelike for A 2 > 1. For A 2 < 1 
the light cones of the metric h are contained in the light cones of g and for A 2 > 1, the 
light cones of g are contained in the light cones of h, as one expects. This is consistent 
with A 2 being the squared sound speed, 



zl = A 2 . 



(B18) 



time 

A 



A < 1 (X > 0) 
/ A 2 > 1 (X < 0) 




lightcone of g 



A > 1 (X > 0) 

A 2 

A < l (X < 0) 



^ space 

FIG. 7: The lightcones of g and h for the given conditions on X and A. The scalar field propagates 
inside the lightcones of h. 

• Spacelike gradient (X < 0) 

In this case, the opposite happens. For A 2 > 1 the light cones of the metric h are 
contained in the light cones of g and for A 2 < 1, the light cones of g are contained 
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in the light cones of h. This is also consistent with our observation that in a static, 
spherically symmetric spacetime radial perturbations propagate with squared speed, 

c 2 s = A" 2 . (B19) 

APPENDIX C: NON-SINGULAR SOLUTIONS: DE SITTER SPACETIME 

Our solutions in the weak field limit of the main text are all singular at the origin r = 0. 
Therefore, our haloes appear to have density cusps at their centers. Of course, near the origin 
the weak field conditions, equations (|13|). are violated, so we do not expect our solutions to 
be valid in this regime. It is interesting to ask whether these singularities are an artifact of 
the weak field approximation, or whether they still persist in solutions of the full non-linear 
equations. Here, by singular we mean that in the coordinate system defined by equation 
(JHJ) quantities like the mass of the halo, the pressure, the energy density and the metric 
coefficients diverge. In this section, we show that the only one non-singular solution of the 
Einstein equations coupled to a k-field whose Lagrangian only depends on X is de Sitter 
spacetime. We note in passing that there are other non-singular static, spherically symmetric 
solutions to these equations, but they typically describe wormholes [4jj instead of haloes, 
and hence, cannot be described in the coordinates we are using. 

Let us begin with the full Tolman-Oppenheimer-Volkoff equations (|15j). We are going 
to look for regular solutions, so we can assume that m, p and p admit Taylor expansions 
around r = 0, 

m = mir + m 2 r 2 + • • • (CI) 
P = Po + Pir + • • • (C2) 
P = Po +Pir H , (C3) 

where we have used that in order for the metric to be non-singular, m = 0. Because both 
p and p depend on X, the squared gradient can be similarly expanded 

X = X + X x r + ■ ■ • . (C4) 

Now, we insert these expansions back into the TOV equations ()15|) and do a term by term 
comparison. Substituting equation (jCl|) into equation (j!5a|) and comparing coefficients we 
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find 

mi = m 2 = 0. (C5) 

Also, inserting equation (|Clj) into equation (|15b|) . one sees than dot/dr is at most linear in r. 
Meanwhile, we substitute the expansions ()C2|) and ()C3|) into the third TOV equation ()15cj) 
and compare coefficients, finding to 0(r _1 ) that p +Po = 0. Therefore, using equation (jSJ), 
we find that either 

Io = or ^ =0. (C6) 
dX x Q 

Using the same TOV equation, and armed with our knowledge that da/dr = 0(r), we obtain 
the relation 2p x + 3pi = 0. Now there are two possibilities; if X ^ 0, then by equation 
()C6|) . pi = —L t x(X )Xi = 0, which implies pi = 0. On the other hand, if X = we find 
that 2pi + 3pi = — 5L t x{X )Xi = 5pi has to vanish, which in turn implies the vanishing of 

Pi- 

We can use these results to set initial conditions at a sufficiently small non-vanishing 
r = Ar. Because po+Po = 0, and because pi = pi = 0, we find that p( Ar) +p(Ar) = 0. With 
this initial conditions we can integrate equation ()15cj) from r = Ar onwards immediately, 
to find that p and p are both constants. This is the only non-singular solution. What is the 
nature of the non-singular solution we have found? Since p is constant and positive for all 
r > 0, the spacetime is simply de Sitter. To verify this is indeed true, we integrate equations 
(I15aj) and (j!5b|) in the following. 

Our result p = — p = constant already solves ()15cj) for all r in one fell swoop. We can 
then integrate (jl5aj) to obtain 

Ait 

m=—r 3 p , (C7) 



which once substituted into equation ()15bj) yields 

1 , (. 8ttG 



a = |log(^l-^lr 2 poj, (C8) 
where we have assumed a(0) = 0. Putting everything back together, the metric is then 

dg2 = ^ _ 87rG r 2 p ^ df2 _ ^ _ ^ r 2 po ^ 1 dr 2 _ r 2 dQ 2 (Cg) 

which is simply de Sitter space in static coordinates. 

It has been known for a while that de Sitter solutions arise for timelike field gradients in 
k-field theories. These solutions exist if the Lagrangian has an extremum at a positive value 
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of X, dL/dX = at Xo > 0. However, for such Lagrangians the term proportional to X 
typically has the opposite sign The constraints that the existence of de Sitter solutions for 
spacelike values of the squared gradient imposes on the Lagrangian are drastically different. 
For example, the conventionally-looking Lagrangian L = X + X 2 does not have a point 
where dL/dX = at X > (timelike), but does have a point where dL/dX = at X < 
(spacelike). Hence, de Sitter solutions at spacelike X might offer an alternative interesting 
possibility to implement inflation and cosmic accelerations with "conventional" Lagrangians. 
In this context, it is instructive to look at our spacelike solution in the coordinate system 
where de Sitter has the form of a flat FRW universe, 



and it can be easily shown that such a (p leads to a constant X = X < 0. The scalar field 
is inhomogeneous, and becomes ill-defined at the horizon, re Ht = H" 1 . 
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